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Abstract 

We propose a new way to construct a multicast coding scheme for linear deterministic relay 
networks. Our construction can be regarded as a generalization of the well-known multicast 
network coding scheme of Jaggi et al. to linear deterministic relay networks and is based on 
the notion of flow for a unicast session that was introduced by the authors in earlier work. 
We present randomized and deterministic polynomial-time versions of our algorithm and show 
that for a network with g destinations, our deterministic algorithm can achieve the capacity in 
[log(.g + 1)] uses of the network. 

1 Introduction 



> 

^ — Computing the capacity and constructing optimal coding schemes for wireless Gaussian networks 

are central open questions and of great importance in network information theory. In a wireless 
network the transmitted signal from a node is broadcasted to all its neighbors and the signal 
• received at a node is the superposition of the signals transmitted by its neighbors and Gaussian 

noise. Broadcasting, interference, and noise are the three main characteristics of a wireless network 
that differentiate it from a wired network and make its analysis much more challenging. Recently 
. . Avestimehr, Diggavi, and Tse [2] proposed an approximation model known as the linear deterministic 

relay network (LDRN) for wireless Gaussian networks that simplifies the three features of wireless 

^ Gaussian networks by instead considering deterministic and linear operations in vector spaces over 

^ finite fields. Avestimehr, Diggavi, and Tse have further shown that for some Gaussian wireless 

networks, the capacity of the wireless network is within an additive constant gap from the capacity 
of the corresponding approximation network and the optimal coding scheme for the approximation 
network can be translated to near optimal coding schemes for the Gaussian wireless network [3]. 

An LDRN is a wireless networking model which can be visualized as a layered directed network 
J\f = {V, E) with set of "nodes" V = IJi^fi where Vi denotes the set of nodes in layer i, and set 
of "edges" E. Let Vi = {vi{\), ■ ■ ■ ,fj(mj)}, where rui denotes the number of nodes in layer i. The 
first layer consists of a single node s = vi(l) called the source node. There are g destination nodes 
denoted by ti = VKi{di),l € {!,••• ,g}, distributed in layers Ki, K2, ■ ■ ■ ,Kg- There is an "edge" 
from every node in Vi to every node in V^+i which corresponds to the transfer matrix between the 
two nodes. Figure [T] is an example of an LDRN with four layers and two destination nodes. 

*The authors are with the Department of Electrical and Computer Engineering, Texas A&M University, College 
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Figure 1: An LDRN with four layers. Here ti = 1^4(1) and ^2 = ^^4(2)- 



During one use of the communication channel between layers i and i + 1, Vi{j) transmits a 
predetermined length vector Xj[j] to the nodes in layer i + 1 and receives a predetermined 

length vector yj+i[fc] given by 

mi 

where Gi[k,j] is a predetermined transfer matrix of the edge {vi{j),Vi^i{k)) £ E. Note that we can 
set Gi[k,j] to be the all-zero matrix if there is no connection from Vi{j) to Vi-^i{k). All vectors and 
matrices are over a fixed finite field F. One can define 
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and the block matrix Gi = [Gi[k,j]] ,1 < k < mj+i, 1 < j < rrii. Then the received vectors at layer 
i + 1 are related to the transmitted vectors at layer i by the following relationship 

yj+i = GjXj. 

The capacity of an LDRN for a single multicast session from source s to the destinations ti, • • • ,tg 
is derived in j2]- Define a cut between the source node s and a destination node tj as a partition of 
nodes V into two sets A and B, with s £ A and tj G B. The capacity of the cut is defined as the 
rank of the transfer matrix from the transmitted vectors of the nodes in A to the received vectors 
of the nodes in 5. [2] shows that the minimum capacity of the cuts between s and tj is the capacity 
of a unicast session between s and tj. Furthermore the multicast capacity of the network between 
source s and destinations ti , • • • ,tg is the minimum of the min-cut capacities between the source 
and each destination. The capacity-achieving scheme in |2] is a random linear coding scheme that 
is asymptotically optimal when the network is used for multiple rounds. 

A few groups of researchers (see, e.g., jH 111 [TTl [12]) have proposed deterministic coding schemes 
for the transmission of a single unicast session over an LDRN which can be constructed in polynomial 
time. Furthermore, they achieve capacity using only one round of the network. These schemes are 
similar to routing schemes in wired networks and have low encoding and decoding complexities at 
the relay nodes. 

In this paper we build upon our work in |1H [T2] to design a simple and low complexity transmis- 
sion scheme for a multicast session over an LDRN. Our scheme will be constructed by progressively 
combining the coding schemes for unicast sessions from the source to each destination. In many 
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ways our scheme is similar to and is a generalization of the scheme in |7] for a multicast session 
in wired networks. We will offer both randomized and deterministic versions of our algorithm and 
show that \log{g + 1)] uses of the network suffice to achieve capacity, which resembles the result 
for wired networks [7J. 

For the case of a single multicast session, there have been multiple recent attempts to devise 
deterministic and efficient algorithms for constructing capacity-achieving coding schemes. In [3], 
Ebrahimi and Fragouli developed an algebraic framework for vector network coding and used this 
framework to devise a multicast transmission scheme over an LDRN. Our scheme has a lower 
complexity of construction and needs fewer uses of the network to achieve capacity. Erez et al. |5] 
offer a different construction by progressing through the network according to a topological order 
and maintaining the linear independence of certain subsets of coding vectors along the processing. 
However, the proposed algorithm does not appear to have a polynomial running time. Kim and 
Medard [9] generalized the algebraic framework of Koetter and Medard [lOJ for classical network 
coding to LDRNs and devised an algebraic algorithm for constructing multicast codes. Again, 
the proposed algorithm does not appear to have a polynomial running time. More recently, [8j 
proposed an algorithm using rotational codes to asymptotically achieve the multicast capacity of 
LDRN networks for a multicast session. Rotational codes have some built-in advantages as they 
are easy to implement at the relay nodes. However, the existence of deterministic polynomial-time 
algorithms for the construction of efficient rotational codes for multicast transmission over an LDRN 
remains unknown. 

We will next review our earlier results on a single unicast session dH [I2] in Section |2] and then 
discuss our coding construction for a multicast session in Section [3] 

2 A single unicast session 

In this section we briefly explain the coding scheme for a single unicast session from |1H 112) . This 
will be the building block of our multicast coding scheme. 

Recall that for each i G {1, ■ ■ ■ , M — 1} the transmitted vector of layer i and the received vector 
of layer i + 1 are related through matrix Gi by yj+i = GjXj. 

For each layer i £ {!,■ ■ ■ , M} label the indices of the vector with the elements of a set Pi and 
label the indices of the vector Xj with the elements of a set Qi . We choose all sets Pi and Qi to be 
disjoint for different values of i. For any A C Pj, let yi{A) denote the subvector of corresponding 
to indices with labels from set A. Similarly, for any B C Qj, let Xi(B) denote the subvector of Xj 
associated with indices with labels from set B. Next partition each set Pi into subset Pj = U™!-^Pj[j] 
and Qi into subsets Qi = ^^iQi[j] such that Pi[j] is the subset of indices of y^ that belong to the 
subvector yi[j] and Qi[j] is the subset of indices of Xj that belong to the subvector Xj[j]. Therefore 
we have yi[j] = yi{Pi[j]) and Xi[j] = Xi{Qi\j]) for any j G {1, • • • , m,} . For any i € {1, • • • , M - 1} 
we will use the sets Pj+i and Qi to label the rows and the columns of the matrix Gi such that for 
each p € P^+i the row of Gi corresponding to the element yj+i(p) is labeled with p and for each 
q € Qi the column of Gi corresponding to the element Xj(g) is labeled with q. For p G Pi+i and 
q G Qi let Gi{p, q) denote the element in row p and column q of matrix Gi. For A C Pj_|_i and B <Z Qi 
let Gi{A,B) denote the submatrix of Gi consisting of the rows in A and the columns in B. Our 
labeling implies that Gi{Pi+i[k], Qi[j]) = Gi[k,j] for any j E {!,••• ,mi} and A; E {1, • • • , mj+i} . 

If node s holds a column vector message w G F^^^ and we are looking at a linear coding scheme, 
then at each layer i G {1, • ' ' ) ; each element of vectors x, and yj will be a linear transformation 
of the vector w. We represent the "global coding vector" (see [7]) for the element Xi{q),q G Qi, with 
row vector Xi{q) € F"*^^^ such that Xj(g) = Xi{q)w and the global coding vector for the element 
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Figure 2: An example of a rate-3 flow from the source node s to the destination node ti. Here the 
matched elements of flow are connected together through dashed lines. 



yi(p))P £ Pi^ with row vector G F-*^^-^ such that yj(p) = yi{p)w. For subsets B (1 Qi and 

j4 C Pj we use the notation Xi{B) and to respectively denote the matrices that are formed 

by the vectors Xi{q) and yi{p) for q £ B and p £ A. Therefore we have Xj(i?) = Xi{B)w and 
yi(A) = 2/i(A)w. 

Suppose that the network supports a rate-i? unicast connection between source node s and a 
destination node t = vxid) for K < M and d G {1, • • • , mi<-} . The main result of [111 [T2] can be 
summarized in the following theorem: 

Theorem 1. For each 1 < i < K and for each 1 < j,k < rui there exist subsets Qi[j] C Qi[j] and 
Pi [k] ^ Pi [k] such that the following hold 

1- \Pi[j\\ = mj]\forie{lr-- ,/^},iG {!,••• 

2- ET=i \P^[j]\ = ET=i =R,foriG{l,--- ,K-1}, 
3. \PK[d]\ = R and \PK[k]\ = for k ^ d, 

4- G'j((J^|^ Pj+i[/c], [jY=i Qi[j]) nonsingular matrix for i £ {1, • • • ,K — 1} . 

Furthermore such subsets can be found by an algorithm that runs in a time that is polynomial in the 
size of the network J\f. 

We call the subsets Qi[j] C Qi[j] and Pi[k] C Pi[k] for i € {1, • • • , M} and j, k G {1, • • • , mj} a 
flow of rate R in the LDRN from the source node s to the destination node t. 

The four properties of a flow in Theorem [T] depend on Gi, Gm-i and do not depend on the 
specific choice of the set Pi[l] among all subsets of Pi[l] with size R. Therefore, if there exists a 
rate--R flow, we can set Pi [1] to be any subset of Pi [1] of size R. 

Notice that the existence of a flow of rate R implies the following simple and low complexity 
coding scheme of rate R from the source s to the destination t: To send message w G F^^^, 
source node s = vi{l) sets yi(Pi[l]) = w and yi(Pi[l] \ = 0. Next, any node Vi{j),i G 

{1, • • • , M} , j G {1, • • • , rrii} , in the network forms the vector Xj[j] by setting 

^^mj])=y^iP^[J])■ 

We say that element p G Pi[j] is "matched" with element q G Qi[j] when Xj(g) is set to yi{p) through 
the preceding equation (see Figure gfor an example of a flow). We further let Xi{Qi[j] \ Qi[j]) = 0. 
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From the properties of flow it follows that at the destination t = vxid) 

^K{QK[d]) = GK-i{pK[d], U Qk-Aj]) • • • GadJ ^[fe], IJ Q2[j])Gi{[J Qi[l])w. 

j=l k=l 3=1 k=l 

Since each matrix Gj(|J^^^ U^i is nonsingular, node t can recover vector w from 

the received vector Xi^((5x[c^]) through a linear transformation. 

3 A coding scheme for a multicast session 

Assume that there are g destination nodes ti, • • • ,tg m the network and the min-cut capacity from 
the source node s to each destination is at least R. We are interested in a multicast coding scheme in 
which all destinations can simultaneously receive the message w G F^^^ of the source. Our scheme 
will be designed by combining the flows of rate R from the source to each destination. 

Suppose that ti = vxiidi) for / € {!,••• ,g}- From Theorem [T] for each ti,l G {!,••• ,g}, there 
exists a flow with subsets P/[A;] C Pi[k] and Q[[j] ^ Qi[j] for I < i < Ki and 1 < j,k < rrii such 
that: 

1- \Phj]\ = \Qi[j]\ for {!,••• ,Ki},j €{!,■■■ 

2. EfJi l^^'bll = ET=i \QM = ^> for i e {!,■■■ ,Ki-l}, 

3. \Pi,^[di]\ = R and = for k / di, 

4. Gi(Ur=t' Pl+M U7=i Qiij]) is a nonsingular matrix for i G {1, ■ ■ ■ , Ki - 1} . 

Since G {I,-- - can be any subset of Pi[l] of size R, we set all subsets G 

{1, • • • , g'} , to be the same subset of 

Our design criterion for a multicast coding scheme is that for each destination ti,l G {!,••• ,g}, 
at each layer i £ {1, ■ ■ ■ , Ki} , the global coding vectors correponding to the elements of the vectors 
yj(P/[j]) for j G {1, • • • ,mj} must be linearly independent vectors and hence the length-i? vector 




can uniquely determine the message vector w. In other words we require for each destination ti and 
each layer i G {1, • • • , Ki} : 

• Condition (*): the matrix (^Uj^i Pi[j]j ™ust be nonsingular. 

The destination node ti = VKi{di) will receive the length-i? vector yKi{PK [di]) = ^^'^(PL 
Since y KiiPR^di]) is a nonsingular matrix, ti will be able to decode message w. 

Notice that at each node Vi{j) for i G {2, • • • , M} we only have control over the design of the cod- 
ing vectors Xi{q) for q G Qi[j] which can be a linear function of the coding vectors {yi{p) : p G -Pi[j]} • 
The coding vectors y^^p) for p G Pi[j] are determined from the coding vectors of the previous layer 
and matrix In our design we will assign coding vectors layer by layer, starting from the first 

layer. At each layer i we fix an arbitrary order on the elements of the set Qi and assign the coding 
vectors Xi{q) in this order. 
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Initialization: We start from the first layer. Since Pl[l] is the same subset for every I E 
{1, • • • ,g} we set yi{P[[l]) = IrxR, i-e., the R x R identity matrix, and set \ PI[1]) = 

for every I S {1, • • • ,g} . In other words we set yi(P|[l]) = w and yi(Pi[l] \ -Pi[l]) = for every 
/ G {1, • • • ,g}. Therefore condition (*) will be satisfied for all destinations in the first layer. 

Inductive Step: We continue our coding construction inductively. Suppose that the condition 
(*) holds for layer i and for all destinations ti = VKi{di) with Ki > i. Next we will design the coding 
vectors Xi(q) for q G Qi one by one and in the order of the elements of Qi so that at the end the 
condition (*) holds for layer i + 1 and all destinations ti = VKi{di) with Ki > i + \. 

At this step of the algorithm for each destination ti with i^; > z + 1 we maintain two matrices. 
One is the matrix Ai which is initially 

(rrii 

and is updated throughout the algorithm. The other matrix is 

Fi = Gi{[j pUMi)^ 

k=l 

where initially Q[ = Uj^iQibl updated throughout the algorithm. Throughout the 

algorithm we maintain the invariance that the product FiAi is a nonsingular matrix for every 
destination ti with Ki >i + \. We will also verify that after all of the elements of Qi are processed, 
for every destination ti with Ki > i + 1 we will have FiAi = yi^i ^Uj=i ^ ^i+i ' "^^lich is sufficient 
for condition (*) to hold at layer i + 

Ai is initially invertible since condition (*) holds for layer i. Matrix Fi is also initially nonsingular 
by the definition of a flow to destination ti given in Theorem[T] Therefore the product FiAi is initially 
nonsingular. Next we will explain the design of the coding vector Xi{q) for q € Qi and describe the 
updating process of Fi and Ai for every destination ti with Ki > i + 1. We consider two cases: 

1. If q is part of the flow for destination ti, i.e., q G Qi[j] for some j G {1, • • • , rrii} , then update 
matrix Ai by replacing row yi{pi) with Xi(q), which we will later explain how to design. Here 
Pi ^ is the unique element that is matched with q G Qi[j] in the flow for destination ti. 
There is no change needed for matrix Fi. 

2. If g is not part of the flow for destination ti, then update Ai adding a new row Xi{q) to it and 
insert a column into Fi so that the set of column indices grows from Q^ to Q'^L) {q}. In this 
step we place Xi{q) in the row of Ai counting from the top which is the same as the position 
of the new column Gi{[J^^^ Pl_^i[k], {q}) in the updated Fi counting from the left. 

When we have gone through all of the elements of Qi, matrix Fi would be ^^(IJ^^^ Pi^i[k],Qi) 
and matrix Ai would be the matrix Xi{Qi). Therefore we have 

FiAi = Gi{ U Pl^[klQi)x,{Q.i) = U Pl^^\j] 

k=i \ j=i 

where the second equation holds since Gi is the transfer matrix from Xi{Qi) = Xi to yij^i- This 
equation guarantees that ^Uj^i^ ^/+ib]) nonsingular, as desired. 

Next we analyze each case and find the condition that Xi{q) needs to satisfy in order for FiAi 
to remain nonsingular: 
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Analysis of Case 1 

Without loss of generality suppose that Xi(q) is the first row of Ai and that matrix Ai after the 
update is of the form 



A, 



Therefore Ai before the update is of the form 
nonsingular. We write 



A' 



ViiPi) 
A' 



We require that the matrix FiAi be 



Fi=[a Fl] 

where ce € F^^^ is the first column of Fi. Using standard matrix calculus we can write 



FiAi =[a F[ 



A' 



cxXi{q) + FlA'i. 



ViiPl) 
A' 



, which is the matrix FiAi resulting from the previous step 



Let us define H = ^ a F[ 

and is nonsingular by the inductive assumption. We can write 

F[A\ = H-ay,{pi) 

and therefore 

FiAi = H + a{xi{q)-yi{pi)). 

For the moment suppose that FiAi is singular. This means that there exist a non-zero column 
vector P e F^''^ with FiAiP = 0. This implies that 

Hp + a{xi{q)-yi{pi))(3 = 0. (1) 

Since H is a nonsingular matrix, there is a vector 'Ji such that o: = H^i. Then ([T| can be rewritten 
as 

Hp + Hjiixiiq) - yi{pi))P = H{P + -fi{xi{q) - yi{pi))P) = 0. 
Since H is nonsingular, the identity holds if and only if 

P + liixi{q)-yi{pi))p = 0. 

If we premultiply the vectors from both sides of the preceding vector equation by {xi{q) — yi{pi)), 
we find that 

{xiiq)-yi{pi))P + {xi{q)-yi{pi))'yi{xi{q)-yi{pi))P = {l + {xi{q)-yi{pi))li){xi{q)-yi{pi))P = 0. 

The expression above is product of two numbers (1 + {xi{q) — yi{pi))^i) and {xi{q) — yi{pi))P- We 
argue that {xi{q) — y^{pi))P is not zero. Observe that if this number was zero, then equation ([T]) 
and the nonsingularity of H would imply that HP and /3 are both zero vectors, contradicting our 
assumption that /3 is a non-zero vector. Therefore 

1 + {Xi{q) - y,{pi))ji = 0. 

This argument implies that for FiAi to be nonsingular it is sufficient to have the following inequality: 

1 + ixiiq) - yM))ll + 0. (2) 
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Analysis of Case 2 

The analysis is very similar to Case 1. Without loss of generality assume that the new row is added 
to the bottom of Ai and the new column is added to the right of Fi. After the update Ai is of the 
form 



Ai 



Here A'^ represents matrix Ai before the update. Also matrix Fi after the update is of the form 

where o: G F^^^ is the new column added to F/, which is the matrix Fi before the update. Our 
inductive assumption implies that H = F[A'i is nonsingular. We can write 

FiAi =H + aXi{q). 

FiAi is singular if there exists a non-zero vector /3 such that 

FiAi(3 = H/3 + axi{q)/3 = 0. 

Since H is nonsingular, there exists a vector 7; such that a = H^i. Therefore FAi is singular if 
there exists a non-zero vector (3 for which 

H(3 + H-fiXi{q)(3 = H{(3 + 'fMq)(3) = 0. (3) 

Since H is nonsingular, ([3| implies 

/3 + 7;a;,(g)/3 = 0. 

If we premultiply both sides of the preceding equation by Xi{q) we obtain 

Xi{q)(3 + Xi{q)-fiXi{q)P = Xi{q)(3{l + Xi{q)ji) = 0. 

The previous equality holds if either Xi{q)/3 = or if (1 + Xi{q)^i) = 0. If Xi{q)(3 = then by ^ 
Hf3 = 0, which, together with the invertibility of H, implies that (3 = 0. But /3 7^ by assumption. 
Therefore if FiAi is a singular matrix, we have 

1 + Xi{q)'yi = 0. 

The preceding argument implies that FiAi is nonsingular if 

l + x,{q)-fi^O. (4) 

A randomized algorithm and the existence of a solution 

Let us summarize the analysis up to this point. The coding vectors Xi{q),q E Qi, can be assigned 
in a way that meet our requirements if 

r= n {l + {x,{q)-y,{pi))^i) J] {1 + x,{q)Jl) ^ 0. 

tl-g&Qi[j],j&{^,- ,mi} ti:q(^Ql[j],je{l,- .mj 

In the preceding equation ti is restricted to the destinations for which Ki > i + 1, and the vectors 
'ji and Viipi) are specified in the analyses of Cases 1 and 2. 
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One other constraint is that Xi{q) for q € Qi[j] can only be a linear combination of the vectors 
iViip) '■ P £ -P«b]} • Let us choose each Xi{q) to be a random linear combination of the elements of 
iViip) '■ P £ -Pib]} where the coefficient of each Uiip) is randomly and independently chosen from 
the uniform distribution over the field F. For each destination ti with Ki > i + 1 define as the 
event that the corresponding term in the product above is zero. Then we have 

Pr(r = 0)=Pr( \/ 0,) < ^ Ft{(Pi). 

ti:Ki>i+l ti:Ki>i+l 

Suppose that q £ Qi[j] and Xi[q) = X^pgp.j^] OpViip)- Now consider a destination ti with Ki > i + 1. 
If g G QilJ] and pi G is matched with q, we need to have 1 + {xi{q) — yi{pi))fi / 0. There 

exist coo e¥,ojp e¥,p e Pilj], which are determined by and -Ji and satisfy 

1 + {xi{q) - yi{pi))'yi = ivo + ^ uipOp 

pePilj] 

There are two cases to consider. First, if = for all p £ Pi[j], then loq + ^p^p.^j ^pOp = wo is a 
constant independent of G Pi[j\- Furthermore by setting 9p^ = 1 and 0p = for p G Pi\j\ and 
p ^ pi so that Xi{q) = yi{pi), we find that 

Wo = 1 + {Xi{q) - yi{pi))li = 1. 

Therefore in this case Pr((^/) = 0. Next if there exists some p G Pi[j] for which ujp ^ {) then 
ujQ + ^pgp.[j] Wp0p depends on 9p,p G Pi[j]- Since 6p,p G -Pj[i], are uniformly distributed random 
variables over F, wo + X]pePjy] ^p^p is likewise uniformly distributed over F. In this case Pr(0;) = p^. 

Next suppose that q ^ Q\[j]- Here we need to have 1 + Xi{q)^i ^ 0. Following the preceding 
argument, there exist G F, Wp G F,p G Pi[j], which are determined by and 7^ and satisfy 

1 + Xi{q)'^i = Wo + ^ Wp6'p. 

If Wp = for all p G Pi[j]-, then ojq + X^pgp.j^j oJpOp = wq is a constant independent of 9p,p G Pi[j]. 
By setting 6'p = for allp G P^b'] so that Xi{q) = 0, we obtain ujq = l + a;i(g)7; = 1; Pr(^;) = 0. 
If there is some p G Pi[j] for which ojp 7^ 0, then an analogous argument to our earlier one implies 
that Pr(^i) = 

As a result, for any destination ti with if; > i + 1, we have Pr((;/!>;) < |^. Therefore 

Pr(T = 0)< ^^'i'l^^WV 

ti:Ki>i+l I I 

Since we are interested in the event that r 7^ 0, we have 

Pr(r7^0)>l-^. 

Therefore if |F| > g, then Pr(T 7^ 0) > and there is at least one valid solution for Xi{q). This also 
yields a randomized algorithm with probability of success of at least 1 — |^. If we take the size of 

the field to be |F| > 2g then the probability of success will be at least 1 — ^ = ^ for each q G Qi. 
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A deterministic polynomial time algorithm 

We next explain a deterministic algorithm with polynomial running time for the finding vectors 
{q)tQ 6 Qi [j] ■ For each q €z Qi we seek a vector u = Xi (q) which is a linear combination of the 
vectors in {yi{p) : p G such that for any destination ti with Ki > i + 1, if q £ Qlij] and 

Pi £ Piij] is matched with q, then 1 + (u — yi{pi))^i / 0, and if g ^ Q\[j] then 1 + w^i / 0. 
Define the subset of indices of destinations W as 



W 



{/ G {I,-- - ,5} : ET/ > i + G for some j G {!,••• , "^0 , y»(M)7i / o} . 



We can write the conditions that u needs to satisfy as 1 + {u — y^{pi))'yi ^ for / G W and 
1 + w^i 7^ for / ^ and Ki > i + 1. Next we use [3 Lemma 8] : 

Lemma 2. Letn < |F|. Letai, ■■■ ,an e ¥^^^ andbi, ■■■ ,bn e F^^^ withaibi 7^ 0, z G {1, • • • , n} . 
There exists a linear combination c 0/ ai, • • • , a„ such that cbi 7^ 0, « G {1, • ' ' > • ^lic/i a vector 
c can be found in time 0{v?R). 

If we are given the set of vectors 7;, then it takes 0{gR) steps to form the set W. Then by applying 
the preceding lemma, if g < \¥\, we can find a vector w G F^^^ such that w is a linear combination 
of the vectors in {yi{pi) : / G W} and for every / G W^, we have that w^i 7^ 0. Furthermore vector 
w can be found in time 0{g^R). By adding the time 0{gR) needed to produce set W, we need a 
total time of 0{g^R + gR) = 0{g'^R) to find vector w. Next we let u = aw for some a G F. We 
show that an appropriate value of a exists such that u satisfies all of the constraints. 

For I & W, we need to have 1 + (aw — yi{pi))ji ^ 0- Therefore 

yMli^. (5) 

For I ^ W and Ki > i + 1 we need to have 1 + crw^i 7^ 0. If W'fi = then this condition is fulfiled 
for all values of a. Otherwise we need to have 

a 7^ — . (6) 

There are at most g constraints of the form (|5]) and ([6| on a. Therefore if the size of field F is 
greater than the number of destinations g, this deterministic approach will find at least one a 
that is not in the discriminating set by at most considering g elements of F. Therefore the total 
complexity of finding an appropriate value of a is 0{g) and the total complexity of finding vector 
u is 0{g + g^R) = 0{g^R). 

To find the overall complexity of finding the vector Xi{q), we need to evaluate the complexity 
of finding vector 7^ for every / G {!,••• ,g} with Ki > z + 1. From the analysis of Cases 1 and 2, 
7; = H^^a, where matrix H is FiAi from the previous step of the algorithm. Since matrix Fi has 
size R X L and matrix Ai has size L x R for some R < L < \Qi\^ computing H needs 0{R\Qi\) 
operations. Evaluating also needs 0{R?) steps and so there are a total of 0{R\Qi\ + ii^) 

operations for evaluating 7^. Since there are at most g different / G {1, • " " i^} with Ki > z + 1, 
we will have 0{gR\Qi\ + gR^) as the total complexity of evaluating different values of 7; for any 
specific q G Qi- Therefore the total complexity of evaluating Xi{q) will be 0{gR\Qi\ + gR? + g^R). 

Let us assume that the number of nodes nii at each layer i G {1, • • • , M} is at most m. Further- 
more assume that the size of transmitted and received signals at each node is at most r. Therefore 
the total complexity of evaluating each Xi(q) will be 0{gRmr + gR^ + g^R). Since there are at most 
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mMr different Xi{q) to be evaluated, if we assume that the unicast flows from source to each desti- 
nation is provided, the total complexity of our algorithm is 0{gRm?' Mr'^ + gB?mMr + RmMr) = 
0{gRmMr{mr + + g)). 

The complexity of computing a unicast flow to a destination by the algorithm given in ^ is 
0{M{mr)^\ogmr). Since we have g destinations, the total complexity of computing the unicast 
flows will be 0{gM(mr)^ log mr). If we add this running time to the running time of our algorithm, 
the total running time will be 0{g'mMr{mrR + R^ +gR+{mr)'^log'mr)). We can compare it to the 
running time of the algoithm given in 0] which is 0{g{r'^mM+R)^ log{r'^mM+R)+r'^mM(r'^mM+ 
Rf + {gloggRMf) and see that our algorithm is considerably faster. 

Number of network uses to achieve capacity 

We have shown that it is sufficient for the size of the field of operation F of the LDRN to be greater 
than g to guarantee the existence of a multicast coding solution. In general however, the network 
operates over some fixed field which is usually Fp for some prime number p. In order to achieve a 
greater field size, we will use multiple rounds of the network. Here we will argue that if we use the 
network for k rounds, it is equivalent to an LDRN with field of operation F = F^. This implies that 
in order to have a field size at least 5 + 1, it is sufficient to use the network for k = \^\ogp{g +1)] 
rounds. This is an improvement over the number of rounds that is needed for the algorithm of |4j 
which is /c ~ \ogp{gi).ogpg — 1)RM). 

Suppose that the network is used for k rounds and we use the superscript t € {0, • • • , A; — 1} to 
denote the time index that a vector is received or sent. For each z S {1, • • • , M — 1} we have 

y*+i = G,x*, t E {0, . . . , - 1} 

Observe we can use a dummy variable D as the unit delay operator and represent the preceding k 
equations as a single equation 

fc-i fc-i 

t=0 i=0 

Next, notice that X^f=o ^^+1-^* St=o -^i-^* regarded as new vectors in the extension 

field F^ and we can assume that the network is operating in the extension field Fp. Since the 
transfer matrix between the layers i and z + 1 is still Gi and has not changed in the new field, the 
existence of the unicast flow over the original field implies the existence of flow over the extended 
field. Therefore our analysis is valid over any field F^. 

References 

[1] A. Amaudruz and C. Fragouli, "Combinatorial algorithms for wireless information flow," in 
Proceedings 2009 ACM-SIAM Symposium on Discrete Algorithms (SODA 2009), New York, 
January 2009. 

[2] A. S. Avestimehr, S. N. Diggavi, and D. N. C. Tse, "Wireless network information flow: A 
deterministic approach," IEEE Transactions on Information Theory, to appear. 

[3] A. S. Avestimehr, S. N. Diggavi, and D. N. C. Tse, "Approximate capacity of Gaussian re- 
lay networks," in Proceedings of 2008 IEEE International Symposium on Information Theory 
(ISIT), 2008, pp. 474 -478. 



11 



[4] J. Ebrahimi, C. Fragouli, "Vector network coding algorithms", Proceedings of IEEE Symposium 
on Information Theory (ISIT2010), Austin, June 2010. 



[5] E. Erez, Y. Xu and E. M. Yeh, "Coding for the deterministic network model," Information 
Theory and Applications Workshop (ITA 2010), San Diego, Feb. 2010. 

[6] M. Goemans, S. Iwata, and R. Zenklusen, "An algorithmic framework for wireless information 
flow," in ^7t/i Annual Allerton Conference on Communication, Control, and Computing , 2009, 
pp. 294-300. 

[7] S. Jaggi, P. Sanders, P. A. Chou, M. Effros, S. Egner, K. Jain, and L. Tolhuizen, "Polynomial 
Time Algorithms for Multicast Network Code Construction," IEEE Transactions on Informa- 
tion Theory, vol. 51, no. 6, pp. 1973-1982, 2005. 

[8] M. Khojastepour and A. Keshavarz-Haddad, "Rotational coding achieves multicast capacity 
of deterministic wireless networks," in Annual Allerton Conference on Communication, 
Control, and Computing. 2009, pp. 1330-1337. 

[9] M. Kim and M. Medard, "Algebraic network coding approach to deterministic wireless relay 
networks," in Proceedings Allerton Conference on Communications, Control, and Computing, 
Monticello, IL, USA, Sept. 2010. 

[10] R. Koctter, M. Medard, "An algebraic approach to network coding," IEEE/ACM Transactions 
on Networking, vol. 11, no. 5, pp. 782-795, Oct. 2003. 

[11] S. M. S. Tabatabaci Yazdi and S. A. Savari, "A max-flow/min-cut algorithm for a class of 
wireless networks," in ACM-SIAM Symposium, on Discrete Algorithms (SODA 2010), Austin, 
Texas, Jan 2010, arxiv.org c-print 0904.2401, April 2009. 

[12] S. M. S. Tabatabaei Yazdi and S. A. Savari, "A combinatorial study of linear deterministic 
relay networks," in Proceedings 2010 Information Theory Workshop (ITW), Cairo, Egypt, Jan 
2010, pp. 1-5. 



12 



